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Based on Takayama-Lin-Liu-Maki model, analytical expressions for the third-harmonic gen- 
eration, DC Kerr effect, DC-induced second harmonic optical Kerr effect, optical Kerr effect or 
intensity-dependent index of refraction and DC-electric-field-induced optical rectification are de- 
rived under the static current-current (Jo Jo) correlation for one-dimensional infinite chains. The 
results of hyperpolarizabilities under Jo Jo correlation are then compared with those obtained using 
the dipole-dipole (DD) correlation. The comparison shows that the conventional Jo Jo correlation, 
?H ' albeit quite successful for the linear case, is incorrect for studying the nonlinear optical properties 

of periodic systems. 



PACS numbers: 78.66.Qn, 42.65.An, 72. 20. Dp, 78.20.Bh 



' I. INTRODUCTION 

c ' 

The different gauge approaches (p • A and E • r) have been adopted in the theoretical studies of both linear and 
nonlinear optical (NLO) properties for many materials [1-3]. For the current-current ( J J) correlation (i.e., the p • A 
gauge), most researchers tend to interpret the current operator J as the static current current Jo [2-6]. For the 
linear transport theory, though the real part of J Jo correlation causes the zero frequency divergence(ZFD), the 
convergent optical properties such as the linear susceptibility, the absorption coefficient, the linear conductivity, etc 
' could be obtained by using the imaginary part of Jo Jo correlation alone, then by applying the Kramers-Kronig (KK) 
relations on the imaginary part of Jo Jo correlation or including the diagmagnetic term [5]. Hence, the static current- 
current (Jo Jo) correlation is widely adopted in the linear transport theory [5-9] and the ZFD is often considered as a 
harmless technical nuisance and tacitly ignored by most researchers. However, for the nonlinear case, Jo Jo correlation 
[9] encounters serious difficulties and the analytical results for nonlinear optical properties do not even converge. 

Among the polymer studies, theoretical calculations of both linear [6-9] and nonlinear optical properties [4,10-19] 
have been carried out based on the different gauges for the simplest 7r-conjugated polymers such as polyacetylene (PA). 
For polyacetylene, some simple periodic, single electron, and tight-binding approximation models such as Su-Shricffcr- 
Heeger (SSH) [20] and Takayama-Lin-Liu-Maki (TLM) [21] have been established to interpret the experimental results 
[22]. But in both linear and nonlinear calculations of the optical properties under the above models, as we pointed 
out recently [9], there are some discrepancies between the conventional treatments using different gauges. Specifically, 
q . if using the same set of wavefunctions but ignoring the phase difference between both gauges and meanwhile applying 
the static current in p ■ A gauge, we cannot guarantee the equivalence between the two gauges, even though the final 
results look quite similar to each other qualitatively. By the example calculation of linear susceptibility under SSH 
model for one-dimensional infinite chains, we strictly proved the nonequivalence between two gauges and ZFD could 
be resolved by considering the gauge factor [9] . Since one needs to apply fairly complicated techniques to resolve ZFD 
in JoJo correlation and preserve the equivalence between two gauges, we prefer the dipole-dipolc(Z3Z3) correlation 
(i.e., the E • r gauge) for nonlinear optical calculations for the polymers. 

On one hand, the DD correlation is derived by assuming a scalar potential E • r as perturbation, giving rise to the 
external electric field E. On the other hand, Jo Jo correlation is obtained by treating the time-dependent uniform 
vector potential A as perturbation. As long as one uses periodic boundary conditions, the scalar potential shows 
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saw-shaped behavior and therefore the resulting electric field is not uniform in the real space, while Jo • A is uniform 
in real space. From this point of view, the Jo Jo correlation seems more appropriate than the DD correlation. Thus 
it is our interest to study some cases which avoid the ZFD difficulties in the J Jo correlation and reveal the pitfalls 
of the Jo Jo correlation via a detailed comparison between DD and Jo Jo correlations. 

Fortunately, the TLM model is one typical case that avoids the ZFD problem, although its sibling model - the SSH 
model is not [11]. The static current operator J derived from TLM model could give us the convergent results for 
hyperpolarizabilities when the frequency approaches 0. However, we consider this result as a mere coincidence, since 
the linear susceptibility calculation based on TLM model diverges in the real part of Jo Jo correlation [9]. Nevertheless, 
we could use TLM model as a common ground to do the comparison between DD and Jo Jo correlations. 

In [23] , we have computed the analytical forms of hyperpolarizabilities for infinite chains by DD correlation under 
both SSH and TLM models. In this paper, we first present a brief description of the static current operator J for 
both models and general formulas for hyperpolarizabilities under J Jo correlation in Section II. We then proceed 
to carry out analytical calculations for DC Kerr effect, DC-induced second harmonic generation, optical Kerr effect 
and DC-electric- field- induced optical rectification by Jo Jo correlation under TLM model for infinite chains (Section 
III). A detailed comparison of the results between DD and Jo Jo correlations is followed subsequently (Section IV). 
The comparison shows that though there are some similarities for some features such as resonant peaks and general 
shapes between these two correlations, important and evident differences abound. For instance, while DD correlation 
clearly indicates the nonexistence of the two-photon cusp in the third-harmonic generation(THG) spectrum [12,13], 
such cusp appeared in Jo Jo correlation; and while DD correlation obviously shows the break of overall permutation 
and Klcinman symmetries [23], Jo Jo correlation maintains both the overall permutation [3] and Klcinman symmetries 
[24] for all frequencies. Finally, we present our conclusions in Section V. 



II. THEORY 



A. nonlinear optical susceptibility under current-current correlation 

The nth-order nonlinear optical susceptibility under current-current (J J) correlation is conventionally reduced to 
the static current-current (Jo Jo) correlation and defined as follows [3,4,9,12]: 

vWro-M a n(e)e2 f i xS„(";"i, ■■-,"») 

n 

with n = — V^Wj, the electronic density, m the electron electron mass, e the dielectric constant, I the unit 

i=l 

dyadic, 5 n> i the Kronecker symbol, and 

1 



n 



W 



j ilv ■ ■ ■ dv„ j d '(//,, 



J dvdte- ik - r+int (fj (r, i)J (ri, h) ■ ■ ■ J (r„, t n )), (2.2) 

where V is the total volume, T is the time-ordering operator, and Jo is the static current operator. 

The Feynman diagram of \^ is simply described as one connected circle in the preceding paper (see Fig. 1 in [23]). 



B. static current operator under SSH and TLM models 

The static current operator J could be found by the commutator between the dipole operator and Hamiltonian. 
For both SSH and TLM models, the current operators were derived in many previous works [4,6-9], here we only list 
the final results. 

For the SSH model, under the same notation of the preceding paper [23], the static current operator Jo.ssh is 
defined by the formula 
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Jq,ssh — — i 

Ls 



e 

ft 



*o + (-l) 



A 



[a - 2(-l) l u] 



i+l,sA,s 



(2.3) 



where to is the transfer integral between the nearest-neighbor sites, A is the gap parameter and C\ s (Ci, s ) cre- 
ates (annihilates) an 7r electron at site I with spin s. 

For the TLM model (Eq.(2.2) in [23]), by adopting the notation in Maki [6] and Wu [4]'s work, the static current 
operator Jo.tlm is defined by the formula 



Jo 



TLM 



v F ^(x)a 3 ^{x), 



(2.4) 



where &(x) = (^{(x) , ^\(x)) is the two-component spinor describing the left-going and right-going electrons, vp is 
the Fermi velocity and a are the Pauli matrixes. 

As pointed out in our recent work [9] , detailed calculations show that the above static current operators lead to the 
ZFD in the linear response for both models. However, in the subsequent calculation for x > we show that the static 
current operator J gives the convergent results for the TLM model. This provides us a convenient base to carry out 
the comparison of the analytical results of \^ between DD and Jo Jo correlations. Hence, the following calculations 
are based on the TLM model only. 



III. HYPERPOL ARIZ ABILITIES FOR TLM MODEL UNDER STATIC CURRENT-CURRENT 

CORRELATION 

A. General four-wave-mixing results 

We apply the general definition Eq.(2.1) and Eq.(2.2) to the TLM model and obtain the following expression for 
Xt1m(^ = +ti, 2 + W3);wi,W2,W3) (or x (3) (^i, u)2, w 3 ) for short): 

(3) 



where yf] is defined by the formula 
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Xjoh h 3 3!L 



]T J^TrL 3 G(k,u ) )a 3 G(k,oj-u )1 ) 



£73 G(k, LJ — U>1 — W 2 )cr 3 G(fc, U — W1—W2 — W3) J , (3.2) 

= 2e h3 ® VF ^S{u ll uj2,u 3 ) (3.3) 

k 

with L the chain length, S(u>i, u>2, 0J3) the summation of the permutations for u>i,u>2,oj 3 , and Green's function G 
defined by the formula 



u) — s + v F k<r 3 + Aai/h 
(uj - s) 2 -ul + ie 



G ( fc > S ) = ~ .,2 , .-. • ( 3 - 4 ) 



In Eq. (3.4), Uk is defined by the formula 

= [{v F kf + (A/ft) 2 ] 1 / 2 . (3.5) 
We now introduce the following three new variables: 

c := A/ft, (3.6) 
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,_!S_,/i + <!*V 



(3.7) 



2^ ~~ 2A' 



(3.8) 



Combining Eq. (3.3), Eq. (3.6)-Eq. (3.8) and replacing the summation over k by its continuous limit, we obtain 



(3) / \ e 4 n «p 

^ (W1,W2,W 3 ) = — ^3- 



2e noVj 
irhA 2 



S{u\,W2,^s)dk 

) 



(3.9) 



(c 5(a>i,W2,w 3 )). 



Substituting Eq. (3.9) into Eq. (3.1), we have 



X (3) (wi,W 2 ,W 3 ) 



2 3 7rA 6 -i(zi + Z 2 + Z 3 )Z1Z2Z 3 J 1 y/ X 2 - 1 



= (c 3 5(wi, u 2 , 0J3)), 



where 



~2~K' 



1...3. 



(3.10) 



(3.11) 



Eq.(3.10) is the general formula for four- wave-mixing under Jo Jo correlation. This is the same as defined in Wu's 
work [4]. As for nonlinear optical susceptibilities, there is no non-equilibrium situation involved, the usage of Kelydish 
Green function in Wu's work is not necessary. 

Now Eq.(3.10) is simplified to compute S(ui, u>2, 1^3) term. In this work, for the purpose of comparing nonlinear 
response between different gauges, we only obtain the analytical formats for third harmonic generation(THG), DC Kerr 
effect (DCKcrr), DC-induced second harmonic gcncration(DCSHG), optical Kerr effect (i.e., intensity-dependent index 
of refraction (IDIR)), and DC-electric-field-induced optical rectification (EFIOR). The results under DD correlation 
with or without contribution in the corresponding figures are obtained from the preceding paper [23]. 



B. Third Harmonic Generation \ (— 3oj; w, oj, oj) 

Applying the Residue theorem and then using Maple to simplify S(lu,u),cj) in the Eq.(3.10), we obtain: 

c 2 (4(5c 2 -2lo 2 ) 8(c 2 -lo 2 ) A{c 2 -2u 2 ) 



LU 4 LU k { 3(4w 2 - OJ 2 ) 3(c 



oj 2 ) (4uj 2 -9uj 2 ) 
(5-8z 2 ) 8(1 -4z 2 ) (l-8z 2 ) 



2 4 c 3 z 4 x \ 3(x 2 - z 2 ) 3(x 2 -4z 2 ) (x 2 - 9z 2 ) 
Combining Eq.(3.10) and Eq. (3.12), we obtain 

e 4 n (fe F ) 3 1 



X {3) (oj,oj,oj) 



1152ttA 6 z 8 
where the function f(z) is defined by the formula 



{3(1 - 8z 2 )/(3z) - 8(1 - 4z 2 )/(2z) + (5 - 8z 2 )/(z)} 



m = 



dx 



1 (a; 2 — z 2 )\J x 2 — 1 



arcsin(z) 

zVT^z 2 

cosh _1 (z) ^ v 

Z\/ 'z 2 — 1 2zy/z 



(z 2 < 1), 

(z 2 > 1). 



(3.12) 



(3.13) 



(3.14) 
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As z — > 0, we have 



m, x e 4 n (hv F ) 3 , 4 32 2 128 , 18944 fi ^ 8xx ,„ ^ 

X (3) = (jg + ^ + — * 4 + -gg-* 6 + 0(z 8 )) (3.15) 



Eq.(3.13) is exactly the same as Wu's result [4]. By applying the following conversion between SSH and TLM 
model: 



then defining 



Hvf = 2t n a, (3.16) 



(3) _ 8 e 4 n (2t a) 3 
Xo "45 ' ^- i7J 



and choosing the same parameters as in the previous works [12,13,23], i.e. A = 0.9eV, no = 3.2 x 10 14 cm 2 and 
a = 1.221, we obtain Xo 3) ~ 1.0 x lO" 10 esu. 

The magnitude of third- harmonic generation under J Jo correlation, and that under DD correlation with or without 
intraband contribution are plotted in Fig.l. The theoretical discrepancies of THG under different gauges have been 
noticed by many others' works [12-19]. It has been addressed in all works that the two-photon absorption peak 
observed in the experiments can not be explained by the single-electron models like SSH and TLM models. Meanwhile, 
it was also pointed out that both gauges should reach the exact same results if the calculations have been performed 
correctly. The reason why the discrepancy exists in the different gauges has not been pinpointed in all previous 
calculations. Recently based on the same models for the linear response [9], we strictly proved that the gauge phase 
factor, which was ignored in the previous studies of optical properties, is the cause for the difference. When the gauge 
phase factor is considered, the difference between different gauges could be resolved [9]. 



C. Optical Kerr effect or intensity-dependent index of refraction x' 3 ' (— w; uj, — u), u>) 

Following an almost identical procedure of obtaining S(u),u),w), we have: 

Be 2 (-48w£ + 60^ 4 c 2 + 24^2 _ 35^2^2 _ 3^4 + 



S(uj, —u),uj) 



3 ul{ut-^){Aul-^f 

1 (-12.x 6 + 15.x 4 + 24a; 4 z 2 - 35a; 2 z 2 - 12a; 2 z 4 + 8z 4 ) 



6c 3 x 3 (x 2 — Az 2 ){x 2 — z 2 ) 3 



Following a similar procedure of evaluating x^(0>0,w), we obtain the optical Kerr effect x^ {—uj; uj, —uj, uj) (or 
X^\uj, — u),uj) for short) as follows: 

(3) e 4 n (fa F ) 3 1 f 2 z 2 (4z 2 -l) 8z 6 -12z 4 + 9z 2 -2 f/ 1 



As z — > 0, we have 



m, x e 4 n (ftw F ) 3 , 4 32 2 128 . 3584 fi fiNX , n nnS 

X< 8 >(", -uj,uj) = -^(^ + ^ + ^ + ^ + 0(z*)) (3.20) 



The magnitude of optical Kerr effect (i.e., intensity-dependent index of refraction(IDIR)) under J0J0 correlation, 
and that under DD correlation with or without intraband contribution arc plotted in Fig. 2. 

Eq.(3.19) is exactly the same as Eq.(13) in Wu's work [4]. From Fig. 2, the results from DD and Jo Jo correlations 
all show the cusp z = 1/2. We would like to point out this is merely the van Hove singularity [8] by the singular state 
density in one-dimensional polymer structure [22], not the real resonant peak. Furthermore, the calculation through 
DD correlations by dropping Vfc terms does not exhibit the z = 1/2 cusp, showing that the cusp is related to the 
process of intraband-transition. 
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D. DC Kerr effect x (3) (-w; 0, 0, w) 



To evaluate the DC kerr effect w;0,0,w) (or x^ (0, 0, w) for short), we first evaluate S(wi, W2, W3) = 

S , (zi,z 2 ,Z3) for general (i = 1,2,3), then substitute it into Eq. (3.10) and take the limit z\ — > 0, z 2 — > 0, 
and z 3 — > z. We obtain 

* (3) (0, 0, «) = e4 5 7 6 ( ^ F 6 )3 zi{ J_ 1)3 {3(3 - 8z 2 )/(z) + (16z 6 - 40z 4 + 18z 2 - 9)} (3.21) 
As z — > 0, we have 

r^w„ „ \ e 4 nn(hvp) 3 , 4 16 9 32 1 256 fi s .. ,„ 

X^(0, 0, W ) = -^-(^ + ^ 2 + ^ 4 + ^ + 0(z 8 )) (3.22) 

The magnitude of DC Kerr effect(DCKcrr) under J J correlation, and that under DD correlation with or without 
intraband contribution are plotted in Fig. 3. The figure only shows one resonant peak at z = 1 for all 3 cases. 

E. DC induced second harmonic generation x^ (— 2o>; 0, u, cu) 

Following a similar procedure of evaluating x^(0>0, w), we obtain the DC induced second harmonic generation 
X^ 3 K— 2o>; 0, lo, uj) (or x( 3 )(0,u;, w) for short) as follows: 



,(3) 

As z — > 0, we have 



, e 4 n (hv F ) 3 1 f2-9z 2 + 2 _ g 4\ 2 j 



i"Vw„ x e 4 no(S.w_F) 3 , 4 16 9 32 4 11776 fi Snn ,„„,n 

X (3) (0,^) = ( T , + ^ + + + 0(z*)) (3.24) 

The magnitude of DC induced second harmonic generation(DCSHG) under J Jo correlation, and that under DD 
correlation with or without intraband contribution are plotted in Fig. 4. The Figure clearly shows two resonant peaks 
at z = 1/2 and z — 1. The width of z = 1 peak suggests that the peak will not be so huge under J J correlation 
than DD correlation if the damping is included. 

F. DC-electric-field-induced optical rectification x^ 3 ' (0; w, — w, 0) 

After the calculations, we obtain the same results as those in DC Kerr effect. Klcinman symmetry [24] is preserved 
in this calculation for all regions. This result is different from DD correlation since Jo Jo correlation maintains the 
commuting feature for all operators. Due to the nonequivalence between EFIOR and DCKerr under DD correlation, 
we still plot the magnitude of EFIOR under Jo Jo correlation, and that under DD correlation with or without intraband 
contribution in Fig. 5. 

IV. DISCUSSIONS 
A. Nonequivalence between DD and JqJq correlations 

From the above calculations in Sec. Ill, the nonequivalence of hyperpolarizabilities between DD and Jo Jo correla- 
tions can be found in all results, though there are some similarities in the resonant peak, the shape of the curve, etc. 
To understand the difference between the gauges in the models, we present a possible explanation in our previous 
work [9]. To maintain the self-completeness of this work, we also briefly address the explanation here. 

If the electromagnetic field is applied, the Schodingcr equation is given by: 



i-(p-gA) 2 + F(r)+# 



V(r,i), (4.1) 
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where ip(r,t) is the exact wave function at space position r and time t, m is the particle mass, q is the electrical 
charge, V(r) is the potential, and A and (f> are vector and scalor potential, respectively. Suppose now A and </> undergo 
the following transformation: 



j A' = A + V/(r,t) 

where /(r, t) is arbitrary, and A' and <j>' are new vector and new scalor potentials after the transformation Eq.(4.2). 
Then it can be shown [25] that the form of the Schodinger equation will be exactly the same if the old wave function 
ip makes the following change into the new exact wave function ip': 

ip -» V' = e iF ^il> = T G (r, t)^, (4.3) 

where gauge phase factor F g (r,t) is defined as: 



F fl (r,t) = |/(r,t). (4.4) 



The above Eq.(4.2) and Eq.(4.3) are called gauge transformation (or U(l) transformation [26]). 

By utilizing the long-wavelength approximation [3,5], the electric field E is described as E = Eoe -4 "*. 

If we consider the following initial scalor and vector potentials under E • r gauge: 

A = 0, (/> = -E • r. (4.5) 

After choosing the gauge phase factor F g as 

by Eq.(4.2), we obtain the new vector and new scalor potential under p • A gauge as: 

A' = -,0' = O. (4.7) 

IL) 

The connection between the old and new wave function is determined by Eq.(4.3). 

Under perturbative schemes to study the optical response, conventionally people use the exact same set of unper- 
turbed wave functions ip^(r,t) °f Hamiltonian Hq (when A = and <ft — in Eq.(4.1)) to serve as our expansion 
basis for both E • r and p • A gauges [2,3]. However, we should point out that the wave functions for both E • r and 
p • A gauges (before and after gauge transformation) should also be restricted by the gauge phase factor F g from 
Eq.(4.3), therefore two basis sets for both gauges are not the exact same unperturbated wave functions ^>„(r, t), but 
arc different by the gauge phase factor F g . And the Hamiltonian under two gauges (E • r and p • A) are not necessary 
equivalent if they are treated independently and are isolated from the connection between the wave functions under 
the two gauges. Unfortunately, this crucial point has not been clearly illustrated and obviously missed by previous 
works using perturbation schemes [1-3]. Especially under current-current correlation scheme, the gauge phase factor's 
contribution is obviously ignored and A 2 (t) term is considered of no physical meanings [3] . Thus the current-current 
correlation is conventionally reduced into the JqJo formula such as Eq.(2.1), and the equivalence between current- 
current and dipole-dipole correlations is usually considered as Jo Jo and DD correlations under the exact same basis 
of unperturbed wave functions [2-4]. 

Langhoff, Epstein and Karplus covered the topics of time-dependent perturbative theory [27] and sharply pointed 
out that the time-dependent phase in wave functions is very essential, the improper treatment of time-dependent 
phase will cause secular divergence in time-dependent perturbations. In field theory, it is also well-understood that 
the improper treatment of the phase factor will cause divergence [5]. Since the gauge phase factor Eq.(4.6) is obviously 
time-dependent, neglecting this phase factor will cause the ZFD in the susceptibility computations. 

Generally speaking, the widely- adopted conventional formula under J Jo is incorrect. It ignores both the gauge 
phase factor's influence and diamagnetic term's contribution [9]. For the linear case, we strictly proved that after 
taking the consideration of the diamagnetic term and the gauge phase factor, both DD and J Jo correlations yield 
the exact the same result for both SSH and TLM models. The details could be found in [9] and will not be repeated 
here. But for the nonlinear case as we mentioned in the THG calculations [12,13], the complexity to include the gauge 
phase factor in J J correlation suggested DD correlation may be more suitable for further studies. 
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B. Zero frequency divergence (ZFD) 



In general, Jo Jo correlation leads to the zero frequency divergence in the nonlinear optical studies. The static current 
operator in the TLM model coincidcntally avoids ZFD problem in the nonlinear calculations shown above, which does 
not mean that it is flawless. For example, linear calculation based on J in TLM model leads a ZFD problem [9]. By 
splitting the Jo term into inter- and intra-band currents in the TLM model and performing the nonlinear calculations 
to determine the contributions from two different currents, we find that the hypcrpolarizabilities for both cases have 
ZFD. For the SSH model, the static current operator J [8] leads to the ZFD in nonlinear calculations [11]. If the 
gauge phase factor [9] is properly considered in our calculations, the ZFD problem could be resolved. Therefore, the 
ZFD problem for nonlinear calculations under the conventional schema of Eq.(2.1) and Eq.(2.2) is not just a technical 
nuisance which can be tacitly ignored. 

C. The overall permutation and Kleinman symmetries 

Based on the p • A gauge, the general formulas of Jo Jo correlations [3,4] preserve both the overall permutation [3] 
and Kleinman symmetry [24] of hypcrpolarizabilities in any systems. Without surprise, our calculations of hyper- 
polarizabilities under J Jo correlation preserve both the overall permutation and Kleinman symmetries. However, 
the overwhelming majority of recent experiments on various physical systems generally refute Kleinman symmetry 
[28]. Based on E • r gauge and ID periodic models, we analytically showed the break of overall permutation and 
Kleinman symmetry [23,29]. Therefore, the experimental testing on the overall permutation symmetry in periodic 
systems can also be used as a valid test for the conventional Jo Jo correlation and p • A gauge. Detailed discussions 
of the symmetry break and some suggested experiments could be found in [23,29]. 

V. CONCLUSIONS 

For the infinite chains under TLM model, the analytical solutions of THG, DCKerr, DCSHG, IDIR and EFIOR are 
obtained through J Jo correlation. The results are not equivalent to those under DD correlations [23]. It shows that 
the conventional Jo Jo correlation formula is incorrect for studying the nonlinear optical properties. Considering the 
complexity of including the gauge phase factor and other terms for the current-current correlation, DD correlation 
may be much more suitable in the nonlinear studies. 
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FIG. 1. The magnitude of third-harmonic generation(THG) under JoJo correlation(dashed-line), under DD correlation 
with(real line) or without(dot-dashed) intraband contribution is in unit of 10 _9 esu. z = hw/2A. 
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FIG. 2. The magnitude of optical Kerr effect or intensity-dependent index of refraction(IDIR) under Jo Jo corrcla- 
tion(dashed-line), under DD correlation with(real line) or without(dot-dashed) intraband contribution is in unit of 10 _9 esu. 
z = hu/2A. 
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FIG. 3. The magnitude of DC Kerr effect(DCKerr) under Jo Jo correlation(dashed-line), under DD correlation with(real 
line) or without (dot-dashed) intraband contribution is in unit of 10 esu. z = ftw/2A. 
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FIG. 4. The magnitude of DC induced second harmonic generation(DCSHG) under Jo Jo correlation(dashed-line), under DD 
correlation with(real line) or without(dot-dashed) intraband contribution is in unit of 10 _9 esu. z = hw/2A. 
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FIG. 5. The magnitude of DC-electric-field-induced optical rectification (EFIOR) under Jo Jo correlation(dashed-line), under 
DD correlation with(real line) or without (dot-dashed) intraband contribution is in unit of 10 -9 esu. z = hco/2A. 
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